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1.  General  Description  of  the  Program  Commencing  November  1.  1988, 

and  continuin  to  the  present,  the  principal  investigator,  his  assoc¬ 
iates  and  research  assistants  have  carried  out  a  program  of  research 
in  the  general  area  of  distributed  parameter  systems  modelling,  iden¬ 
tification  and  control.  This  work  has  been  concentrated  in  four  main 
areas : 

1)  Control  of  the  Korteweg  -  de  Vries  and  other  nonlinear  par¬ 
tial  differential  equations  exhibiting  solitary  waves; 

il)  Frequency  domain  analysis  of  distributed  parameter  systems; 
ill)  Energy  dissipation  mechanisms  in  elastic  structures; 

Iv)  Modelling  and  parameter  estimation  in  elastic  structures. 

The  research  program  has  included  mathematical  studies  by  the  princi¬ 
pal  investigator,  contributions  from  visiting  associate  researchers 
and  dissertations  prepared  by  PhD  students  of  the  principal  investi¬ 
gator.  Related  activities  include  experimental  work  in  the  .MIPAC-VT 
laboratory  by  the  principal  investigator  and  visiting  associates, 
presentation  of  results  at  scientific  meetings  and  preparation  and 
publication  of  research  articles.  In  the  sequel  we  first  provide 
brief  descriptions  of  each  area  of  research  and  the  results  obtained 
and  then  summarize  research  publication  activity,  visitors,  scientif¬ 
ic  meetings,  etc. 
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2 .  Brief  Descriptions  of  Research  Areas  and  Results 

i  )  £ont£o_l  lie_K£r_t ewe£  ^  de_V£i£S_and_o^h£r_n£nl^i near_p^r^ial 

dif.f ere^n£lal_e£uatj^ons_exhib^t_in£  ^ol^i^a£y_wav£s  . 

With  the  control  theory  of  linear  partial  differential  equations 
having  now  reached  a  certain  state  of  maturity,  we  have  felt  it  to  be 
Important  to  initiate  studies  for  nonlinear  systems  with  significant 
applications.  As  in  the  case  of  the  early  development  of  the  control 
theory  for  linear  distributed  parameter  systems,  we  have  found  it  de¬ 
sirable  to  focus  on  nonlinear  systems  for  which  a  considerable  body 
of  theory  exists  in  regard  to  the  classical  questions  of  existence, 
uniqueness  and  regularity  before  in  order  to  obtain  definitive  cont¬ 
rol  results.  At  the  present  time  we  are  studying  two  nonlinear  dis¬ 
persive  wave  equations,  the  Korteweg  -  de  Vries  (KdV)  equation  and 
the  related  Bousslnesq  equation  and  a  nonlinear  diffusion  equation, 
the  Benjamin  -  Bona  -  Mahony  (BBM)  equation.  A  number  of  important 
positive  results  have  been  obtained  in  the  KdV  and  BBM  cases.  It  is 
significant  that  the  study  of  these  equations  indicates  questions  for 
related  linear  systems  which  otherwise  might  be  ignored.  Here  we 
will  describe  results  obtained  for  the  KdV  equation;  research  is  in 
preliminary  stages  for  the  other  two. 

Work  carried  out  jointly  with  Professor  B.-Y.  Zhang  of  the  Uni¬ 
versity  of  Cincinnati  and  Prof.  V,  Komornik  of  the  University  of 
Strasbourg,  France,  involves  the  forced  KdV  equation 

.  w  M  =  f  (1.01) 

3t  3x  3x® 

on  the  region  0  <  x  <  2r,  t  >  0,  with  periodic  boundary  conditi- 
tlons 

w^*''  (2:T,t)  =  w^^^  (O.t)  .  k  =  1,2,3 . 

so  that  X  =  27:  is  effectively  identified  with  x  =  0. 

For  the  unforced  equation  (f(x.t)  ■  0)  it  is  well  known  that 

there  is  an  infinite  set  of  conserved  quantities.  Only  the  first 
three  of  these  are  of  Interest  to  us: 
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f  w(x.t)  dx  .  (1.02) 

«  0 

.2t: 

w(x.t)*  dx  ,  (1 .03) 

‘'O 


(1.04) 


Much  of  our  control  analysis  involves  examination  of  the  behavior  of 
these  quantities  when  the  forcing  term.  f(x,t),  is  not  identically 
zero.  For  suitably  smooth  solutions  of  (1.01)  we  have 


2i: 


,2n 


fl  3 

-I  w(x,t)  dx  =•  f-wCx.t)— (x.t)-.^-4(x,t)  +  f  (x.t)] 

atjQ  Jq  l  ax  ax*  J 


dx  = 


'  r  *  I* 

[— f-2W(x,t)*  -  ^-~(x,t)]  +  f(x,t)]  dx  »  f(x,t)  dx  . 


2r 

0 


from  which  we  conclude  that  "volume"  is  conserved  if  we  require 

.2tt 


J 
«  0 


f  ( X  .  t )  dx  =  0  . 


(1.05) 


In  fact  we  are  interested  in  a  particular  sort  of  control 
action.  Let  I  =  [a.b]  be  a  subinterval  of  [0.2-].  require 

f(x.t)  «0  ,  x^I  . 
and  determine  f  on  I  by 

f(x.t)  =  -k(w(x,t)  -  [w]j)  .  (1.06) 


where  k  is  a  positive  "gain"  constant  and 
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[w] 


I 


[w]l(t) 


1 

b-a 


j.b 

a 


w( X . t ) 


dx 


is  the  aean  value  of  w  over  I.  It  is  trivial  to  verify  that  (1.05) 
is  satisfied.  This  sort  of  control  can  be  realized,  approximately, 
by  a  distributed  "pumping”  action  over  the  interval  I  which  injects 
"fluid"  into  the  system  at  points  x  where  w(x,t)  <  [wlj^  b]^^^ 

extracts  "fluid"  from  the  system  where  the  opposite  inequality  ob¬ 
tains,  at  a  proportional  rate  determined  by  k  in  each  case. 

The  feedback  control  law  just  described  is  intended  to  stabilize 
the  controlled  system  in  a  certain  sense.  Since  volume  is  conserved, 
for  a  solution  of  (1.01)  with  initial  state 

w(x,0)  =  Wq(x) 


we  have,  for  all  t. 


V(w(- ,t)) 


2r 


w(x,t)  dx  *  V(w^)  = 


w(x,0)  dx  ; 


equivalently,  abbreviating  [wJ^O  2-]  convenience 


Iw](t)  ■  [w^]  . 


We  then  note  that  for  any  w  e  L®t0.2r]  satisfying  this  Identity 


2t: 

w(x , t ) *  dx 

0 


[w] 


dx  = 


w(x ,  t ) 


.2- 


+  2[w] 


^w(x,t)  -  [w]j  dx 


V 

=  0 


n2r 

[w]®  dx 

0 
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^w(x.t)  -  [w]j  dx  + 

0 

Hence,  among  functions  w  for  which  V(w)  =  V(w^),  the  integral  (1.03) 

is  minimized  by  the  constant  function  [w] .  We  may  therefore  hope 
to  cause  w(x,t)  to  approach  [w]  as  t  -»  <»  by  using  a  control  law  de¬ 
signed  so  that  this  Integral  is  non-increasing. 

Since  the  KdV  equation  has  the  property  of  invariance  under 
time  reversal  (replace  x  by  -x  at  the  same  time  as  we  replace  t  by 
-t),  global  controllability  results  for  the  KdV  equation  are  in  our 
grasp  if  we  can  use  the  LaSalle  Invariance  Principle  to  obtain  global 
asymptotic  stability  of  the  zero  equilibrium  state  and  if  we  can 
obtain  local  controllability  results  in  a  neighborhood  of  that  state. 

The  application  of  the  LaSalle  Invariance  Principle  referred  to 
above  requires,  for  its  Implementation,  certain  unique  continuation 
results  for  solutions  of  the  unforced  KdV  equation.  We  have,  in  fact, 
shown  that  any  nontrivial  solution  of  the  KdV  equation  lying  in  an 

appropriate  function  space  (l“{-««>,»)  in  this  case)  cannot  have  supp¬ 
ort  restricted  to  two  left  (or  right)  horizontal  lines  in  x-t  space. 
This  result  implies  that  the  solution  u  cannot  have  compact  support 
at  two  different  times  and  also  implies  that  if  the  solution  u  van¬ 
ishes  in  an  open  set  in  x-t  space,  then  it  vanishes  everywhere.  It 
turns  out  that  these  results  provide  the  proof  of  triviality  of  cer¬ 
tain  invariant  sets  for  the  closed  loop  KdV  system  described  above 
which  is  essential  for  completion  of  the  asymptotic  stability  result. 
In  order  to  complete  the  global  asymptotic  stability  part  of  the  work 
it  is  necessary  to  show  that  solutions  of  the  controlled  system  lie 

in  a  compact  subset  of  the  state  space,  L*[0,2r].  This  part  of  the 
work  has  proved  remarkably  recalcitrant  and  requires  our  continuing 
strenuous  efforts  for  successful  completion. 

The  required  local  controllability  results  for  the  nonlinear  eq¬ 
uation  in  a  neighborhood  of  w  ■  0  follow  from  exact  controllability 
results  for  the  linearized  system,  which  is  just  the  third  order  lin¬ 
ear  dispersion  equation 


J 


2r 

[w]®  dx 

0 


3w  ^  3^w  ^  j 
3t  3x* 


*1 
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We  have  conpleted  a  thorough  study  of  the  uniform  exponential  stabil- 
izability  and  exact  controllability  of  this  system  and  a  related 
boundary  control  system  with  positive  results  in  each  instance.  This 
work  is  complete  and  will  appear  as  a  paper  jointly  authored  by  B.- 
y.  Zhang  and  the  principal  investigator.  The  results  on  the  KdV  equ¬ 
ation  are  currently  being  prepared  as  an  article  for  publication 
jointly  authored  by  V.  Komornik,  B.-Y.  Zhang  and  the  principal  inves¬ 
tigator  . 

i  i  )  Fr^queiicy  ^oma  j^n_aiia2.y£i  o^f _d^s_t r^but^d_pa^ra^m^t^r_sis_t ems_^ 

In  this  area  we  have  been  concerned  with  the  development  of  fre¬ 
quency  domain  methods  useful  for  analyzing  hyperbolic  systems  corres¬ 
ponding  to  various  types  of  wave  motion,  including  vibration  of  elas¬ 
tic  systems  of  various  types.  The  transfer  functions  and  related  in- 
out  -  output  operators  corresponding  to  these  types  of  systems  have 
some  disconcerting  properties,  particularly  in  relation  to  converg¬ 
ence  questions  with  respect  to  approximation  via  finite  dimensional 
systems.  In  addition  to  these  approximation  questions  we  have  been 
studying  the  use  of  transfer  function  methods  to  analyze  closed  loop 
systems  arising  out  of  certain  linear  feedback  laws  in  distributed 
parameters  systems,  particularly  where  boundary  feedback  is  concern¬ 
ed,  use  of  transfer  functions  to  determine  the  admissibility  of  var¬ 
ious  input  and  output  mechanisms  and  a  number  of  other  related  mat¬ 
ters.  In  this  process  we  have  developed  a  representation  of  closed 
loop  semigroups  in  terms  of  the  inverse  Laplace  transforms  of  the 
associated  closed  loop  transfer  functions,  providing  what  we  believe 
to  be  a  significantly  new  approach  to  certain  problems  in  this  area. 

Our  work  in  this  area  has  been  greatly  strengthened  by  the  pres¬ 
ence  of  Dr.  G.  Weiss  as  a  visiting  researcher.  He  has  developed  the 
particularly  useful  category  of  regular  linear  systems  which  are,  to 
put  it  briefly,  linear  input-output  systems  whose  response  to  an  im¬ 
pulse  input  occurring  at  t  =  0  has  a  certain  limiting  behavior  as 
t  1  0 .  The  introduction  of  this  class  of  systems  and  the  systematic 
development  of  its  properties  has  clarified  many  previously  obscure 
questions . 

Dr.  Weiss  has  been  investigating  the  connection  between  state 
space  and  transfer  function  representations  of  infinite  dimensional 
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llaear  systems.  The  central  question  Is  to  find  verifiable  condi¬ 
tions  for  a  system  to  be  regular,  since  regular  systems  are  known 
to  have  a  convenient  state  space  representation. 

Dr.  Weiss  has  been  working  with  the  principal  investigator  in 
some  research  on  well  posedness  of  feedback  systems  when  the  feedback 
operator  is  unbounded.  Conditions  have  been  found  for  the  closed 
loop  system  to  be  regular.  Explicit  characterizations  of  the  closed 
loop  generator,  control  operator  and  observation  operator  have  been 
obtained.  More  recent  work  with  the  Principal  Investigator  has  been 
focussed  on  the  question  of  a  frequency  domain  characterization  of 
exact  controllability  and/or  observability  of  infinite  dimensional 
linear  systems,  extending  to  infinite  dimensions  a  criterion  due  to 
M.  L.  Hautus  in  the  finite  dimensional  setting.  This  research  is  now 
complete  and  has  been  submitted  for  publication  in  the  SIAM  Journal 
on  Control  and  Optimization  as  a  Joint  paper  by  6.  Weiss  and  the 
principal  investigator. 

Dr.  Weiss  has  been  cooperating  with  the  principal  investigator's 
former  PhD  student.  Dr.  Scott  Hansen,  on  extension  of  the  Carleson 
measure  criterion  for  admissibility  of  control  operators  to  the  case 
when  the  input  space  is  infinite  dimensional.  They  have  also  invest¬ 
igated  the  connection  between  admissibility  and  the  related  operator 
Lyapounov  equation. 

Working  with  Prof.  Richard  Rebarber  of  the  University  of  Nebras¬ 
ka,  Dr.  Weiss  has  found  some  general  conditions  for  a  feedback  system 
to  be  robustly  stable  with  respect  to  small  delays  in  the  feedback 
loop.  Specifically,  they  have  found  that  to  ensure  robustness,  the 
open  loop  transfer  function  should  decay  in  all  directions  in  the 
right  half  complex  plane  at  a  certain  rate.  This  research  is  also 
complete  and  will  appear  as  a  paper  jointly  authored  by  Dr.  Weiss 
and  Prof.  Rebarber. 


ill)  I.n^r£y_dis^l£a_t l£n_m^chanis.m^  i.n_eias.tic_s t.ruct.ur.e^. 

Our  research  program  in  this  area  has  been  concentrated  on  the 
study  of  internal  dissipation  mechanisms  in  the  context  of  the  Euler- 
Bernoulli  beam  equation 
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(3.01) 


for  which,  in  the  absence  of  additional  dissipative  terms,  the  energy 


■'■S'  -iD'©’ 


(3.02) 


is  conserved.  Difficulties  encountered  in  developing  mathematical 
models  for  energy  dissipation  in  this  system  agreeing  with  dissipa¬ 
tion  rates  measured  in  the  MIPAC-VT  laboratory  led  earlier  to  the 
development  of  an  approach  which,  accurately  or  not,  has  been  named 
the  spatial  hysteresis  damping  model.  This  model  consists,  primarily, 
in  the  Integro-partlal  differential  equation 


2  ^  f  h(x.  f)  ri^(x,t)-^^(f  ,t)]df 
3x  Jq  Lstax  3t3x  J 


(3.03) 


involving  an  appropriate  interraction  kernel  h(x,0*  The  analytical¬ 
ly  derived  damping  rates  associated  with  this  model  show  excellent 
ar'reement  with  laboratory  measurements  over  a  wide  frequency  range, 
as  described  in  reports  filed  earlier. 

Damping  mechanisms  such  as  this  one  are  direct  mechanisms  in 
that  they  Involve  direct  insertion  of  supplementary  dissipation  terms 
into  the  original  conservative  equations  governing  the  elastic  sys¬ 
tem.  Mechanisms  studied  in  our  recently  completed  work  are  indirect 
in  the  sense  that  they  involve  coupling  the  mechanical  equations  gov¬ 
erning  beam  motion  to  related  dissipative  systems  with  additional  dy¬ 
namics,  resulting  in  an  overall  system  in  which  mechanical  energy  is 
dissipated.  We  have  carefully  analyzed  the  Euler-Bernoul 1 i  beam  with 
thermoelastic  damping  and  with  shear  diffusion  damping.  The  most  sig¬ 
nificant  situation  is  that  wherein  both  types  of  damping  are  simul¬ 
taneously  present.  The  paper  embodying  this  research  has  been  comp¬ 
leted  and  has  been  accepted  for  publication  in  the  Quarterly  of  Ap¬ 
plied  Mathematics.  It  will  appear  shortly. 

Indirect  damping  mechanisms  enjoy  the  advantage  that  they  are 
motivated  and  derived  from  simple  physical  considerations,  allowing 
us  to  answer  objections  often  made  concerning  the  spatial  hysteresis 
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■  odel  to  the  effect  that  it  is  ultimately  act  hoc  in  nature.  Because 
of  these  desirable  properties  we  have  thought  it  useful  to  study  such 
mechanisms  from  a  fairly  general  viewpoint.  In  an  article  just  comp¬ 
leted  we  have  introduced  a  general  analytic  and  algebraic  structure 
within  which  many  such  indirect  damping  mechanisms,  including  but  not 
restricted  to  the  thermal  and  shear  diffusion  mechanisms,  may  be  in¬ 
corporated.  The  resulting  framework  is  sufficiently  general  to  Indi¬ 
cate  the  possibility  of  many  other  indirect  damping  mechanisms  than 
those  already  dealt  with  and  precise  enough  so  that  it  permits  us  to 
see  that  thermal  damping,  even  when  longitudinal  heat  conduction  is 
admitted,  is  mathematically,  as  well  as  physically,  quite  distinct 
from  shear  diffusion  damping.  This  article  has  been  submitted  to  the 
Journal  of  Mathematical  Analysis  and  Applications. 

In  this  area  we  have  also  been  concerned  with  the  development  of 
new  laboratory  techniques  which  will  allow  us  to  overcome  experiment¬ 
al  difficulties  long  encountered  in  separating  the  effects  of  inter¬ 
nal  damping  from  those  due  to  interaction  with  the  supporting  appara¬ 
tus  and/or  measuring  devices.  A  new  technique  has  now  been  developed 
allowing  one  to  study  waves  in  the  elastic  structure  as  they  leave 
the  supporting  clamp  structure,  comparing  them  with  waves  returning 
to  that  structure  after  reflection  from  the  free  end  of  the  elastic 
beam  under  study.  This  technique  uses  MIPAC-VT's  laser  vibrometer 
to  eliminate  effects  due  to  measurement  devices,  as  are  common  with 
the  use  of  accelerometers  and/or  stress  gauges,  and  relies  on  the 
mathematically  provable  fact  that  damped  modes  in  an  elastic  struc¬ 
ture  do  not  have  true  nodes  as  developed  theoretically  for  undamped 
structures.  Internal  damping  rates  are  inferred  from  comparison  of 
vibration  amplitudes  at  quasi-nodes  near  the  clamping  mechanism  with 
corresponding  quasi-nodal  amplitudes  near  the  free  end  of  the  struc¬ 
ture,  It  is  these  easily  made  observations  which  allow  experimental 
assessment  of  the  outgoing  and  incoming  wave  amplitudes  as  discussed 
earlier.  A  paper  on  this  experimental  technique,  along  with  some 
experimental  results,  is  currently  in  preparation. 

Iv)  Mo^e^li^n^  ^nd  £a£ame^er  ^s_t inia_t i^n_i n  £l^s_ti£  ^txu£tur£s_;_ 

Here  we  are  referring  to  the  work  described  under  section  c) 
of  the  Renewal  Proposal  to  which  this  appendix  is  attached.  We  will 
not  repeat  the  basic  discription  of  the  research  program  as  developed 
there.  Rather,  we  will  indicate  what  work  has  been  carried  out  up  to 


-10- 


the  present  time  with  the  support  of  AFOSR  89  -  0031.  This  work  has 
been  a  cooperative  undertaking  with  Professor  Luther  White  of  the 
University  of  Oklahoma ,  whose  research  is  also  partially  supported  by 
AFOSR . 

Extensive  use  of  the  research  equipment  in  the  MIPAC-VT  facility 
has  been  made  by  the  principal  investigator,  Assistant  Professor  Rob¬ 
ert  C.  Rogers,  and  Luther  White,  visiting  from  the  Unive.sity  of  Ok¬ 
lahoma,  in  study  of  modal  frequencies,  shapes  and  damping  rates  in  a 
variety  of  elastic  structures.  Major  efforts  have  been  made  to  imp¬ 
rove  laboratory  capabilities  in  this  area  with  the  design  and  instal¬ 
lation  of  equipment  allowing  separate  excitation  of  lateral  and  tor¬ 
sional  modes  in  axially  symmetric  narrow  plates.  This  equipment  was 
first  put  to  use  during  a  spring  1991  visit  by  Luther  White,  with 
experimental  results  far  better  than  any  obtained  previously.  This 
data  is  currently  being  analyzed  for  identification  purposes. 

This  experimental  data  is  complemented  by  our  efforts  to  develop 
usable  narrow  plate  models,  starting  with  the  .Mindlin  -  Timoshenko 
plate  model.  Joint  work  by  the  principal  investigator  and  Prof.  Lut¬ 
her  White  has  now  resulted  in  such  a  model  which  appears  to  offer 
real  advantages  from  the  identification  standpoint.  This  model  is 
a  distributed  model  Insofar  as  the  longitudinal  axis  of  the  narrow 
plate  is  concerned  but  is  finite  dimensional  as  far  as  description 
of  shearing  and  torsional  motions  orthogonal  to  the  longitudinal  axis 
is  concerned.  Mathematically,  the  resulting  model  is  a  linear  sym¬ 
metric  hyperbolic  system;  an  extension  of  the  familiar  Timoshenko 
beam  equations.  It  is  far  simpler  to  work  with  than  the  two  dimen¬ 
sional  plate  equations  and,  in  that  respect,  appears  to  offer  real 
computational  advantages  for  identification  purposes. 


3 .  Research  Associates  and  Assistants  Supported  in  Part  by 
AFOSR  89  -  0031 


In  accordance  with  the  provisions  of  the  subject  grant,  three 
visiting  research  associates  were  supported  during  the  grant  period. 

Dr.  Gunter  Leugering  of  Technische  Hochschule  Darmstadt,  Darm¬ 
stadt,  Germany,  was  supported  in  part  by  grant  funds  during  the  per¬ 
iod  January  1.  1989  through  April  30,  1989  as  a  visiting  Assistant 

Professor  in  the  Mathematics  Department  at  Virginia  Tech.  Dr.  Leuger- 
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ing  assisted  the  principal  investigator  in  development  of  the  shear 
diffusion  damping  model  and  the  related  projection  method.  He  also 
carried  out  research  in  control  theory  of  partial  differential  equa¬ 
tions  with  delay  terms,  such  as  arise,  for  example,  in  connection 
with  elastic  systems  with  viscoelastic  damping. 

Grant  funds  were  used  during  the  period  covered  by  the  grant  to 
provide  partial  support  for  Dr.  George  Weiss,  who  visited  Virginia 
Tech  from  the  Weizmann  Institute,  Rehovoth,  Israel.  It  was  sufficient 
to  underwrite  only  20%  of  Dr.  Weiss  salary  as  the  larger  part  was 
provided  in  the  form  of  a  Weizmann  Fellowship.  Dr.  Weiss  carried  out 
an  intensive  program  of  research  in  the  area  of  frequency  domain  met¬ 
hods  for  infinite  dimensional  linear  systems.  His  appointment  exten¬ 
ded  from  December  15,  1989  through  December  14,  1990. 

Our  research  efforts  in  nonlinear  distributed  parameter  control 
were  greatly  aided  by  the  visit  to  Virginia  Tech,  during  Sem.  I, 

AY  1990-91,  by  Professor  V.  Komornlk,  newly  appointed  Professor  of 
Mathematics  at  the  University  of  Strasbourg,  France.  Professor  Kom- 
ornik,  who  was  partially  supported  by  this  AFOSR  grant,  has  gained 
world  recognition  for  his  exceptionally  insightful  treatment  of  con¬ 
trol  problems  for  wave  and  elastic  systems,  as  well  as  certain  non¬ 
linear  distributed  parameter  systems.  He  worked  with  the  principal 
Investigator  and  the  principal  investigator's  former  PhD  student,  Dr. 
B.-Y.  Zhang  on  global  controllability  of  the  Korteweg  -  de  Vries  eq¬ 
uation. 

The  principal  investigator's  former  PhD.  student,  Scott  Hansen, 
was,  in  part,  supported  by  the  grant  from  November  1,  1988  through 
April  30,  1989.  .Mr.  Hansen  had  also  been  supported  in  part  by  AFOSR 
at  the  University  of  Wisconsin  from  an  earlier  grant  at  that  Univer¬ 
sity.  In  December  of  1989  Mr.  Hansen  received  the  Mathematics  PhD  at 
the  University  of  Wisconsin  -  Madison,  presenting  a  thesis  titled 
Frequency  -  Proportional  Damping  Mechanisms  in  Elastic  Beams.  This 
thesis  dealt  with  mathematical  properties  of  the  spatial  hysteresis 
damping  mechanism  and  related  "bending  rate"  models,  with  particular 
attention  being  given  to  the  properties  of  the  strongly  continuous 
semigroups  associated  with  such  models  and  the  character  of  the  gen¬ 
erating  operator  of  the  semigroup.  Dr.  Hansen  continued  these  re¬ 
searches  at  Virginia  Tech  during  Sem.  II,  1988-89.  He  is  currently 
at  Iowa  State  University. 

The  principal  investigator's  former  PhD  student.  Mr.  B.  -  Y. 
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Zhang,  received  his  PhD  from  the  University  of  Wisconsin.  Madison, 
in  May,  1990,  presenting  a  thesis  titled  Some  Results  for  Nonlinear 
Dispersive  Have  Equations  with  Application  to  Control.  Although 
Mr.  Zhang  was  not  directly  supported  by  grant  funds,  his  thesis  re¬ 
search  nevertheless  parallelled  and  supported  the  principal  inves¬ 
tigator's  work  in  the  control  of  nonlinear  partial  differential  eq¬ 
uations.  Dr.  Zhang  is  now  Assistant  Professor  at  the  University  of 
Cincinnati . 


4 .  Research  Articles  Resulting  from  Work  Supported  in  Part  by 
AFOSR  89  -  0031 


The  following  is  a  list  of  papers  reporting  research  supported 
in  part  by  A*'0SR  89  -  0031.  In  the  case  of  papers  by  the  principal 
investigator  it  includes  some  articles  for  which  the  research  was 
carried  out  under  the  aegis  of  earlier  AFOSR  grants  but  for  which 
processing  continued  under  AFOSR  89  -  0031.  Listed  papers  by  other 
authors  report  only  research  carried  out  under  AFOSR  89  -  0031. 

£a£er,s_b3i^  ^he,  j2.ri.nc.iE.al.  i.nves.ti.gato.r.L 

On  mathematical  models  for  the  elastic  beam  with  frequency 
proportional  damping,  in  "Control  and  Estimation  in  Distributed 
Parameter  Systems,  H.  T.  Banks,  Ed.,  SIAM  Publ.,  Philad.,  1990 

On  the  positive  square  root  of  the  fourth  derivative  operator. 
Quarterly  of  Applied  Mathematics,  46  (1988),  pp .  751  -  773. 

Spectral  and  asymptotic  properties  of  linear  elastic  systems 
with  internal  damping.  Proc.  Conf.  on  Boundary  Stabilization 
and  Control  of  Systems  Governed  by  PDE '  s ,  Clermont-Ferrand,  1988 

(with  T.  P.  Svobodny)  Phase  identification  in  linear  time-perio¬ 
dic  systems.  IEEE  Transactions,  Vol.  AC  34  (1989),  pp .  218-220 

(With  G.  Chen  et  al)  Analysis,  designs  and  behavior  of  dissipat¬ 
ive  Joints  for  coupled  beams.  SIAM  J.  Appl.  Math.  49  (1989),  pp. 
1665-1693 

Computational  study  of  the  Korteweg  -  de  dries  equation  with 
localized  control  action,  in  "Distributed  Parameter  Control 
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Sysfews:  New  Trends  and  Applications" ,  G-  Chen,  E.  B.  Lee, 

W.  Llttman  and  L.  Markus,  Eds.,  Marcel  Dekker,  Inc.,  Sew  York, 
1990,  pp.  195  -  203 

Approximation  of  input-output  operators  for  distributed  parameter 
systems.  Proc.  1990  lEEE-SIAM  Conf.  on  Dec.  &  Cont.,  Honolulu, 
Oecenber  1990. 

A  comparison  of  certain  elastic  dissipation  mechanisms  via  de¬ 
coupling  and  projection  techniques ,  to  appear  in  Quarterly  of 
Applied  Mathematics,  1991 

Neutral  FOE  canonical  representations  of  hyperbolic  systems. 

Jour.  Int.  Eq'ns.  3  (1991),  129  -  166 

A  general  framework  for  elastic  systems  with  indirect  damping, 
submitted  to  J.  Math.  Anal.  Appl . 

(With  G.  Weiss)  A  general  necessary  condition  for  exact  observa¬ 
bility,  submitted  to  SIAM  J.  Cont.  Opt. 

(With  B.-Y.  Zhang)  Controllability  and  stabi lizabi lity  of  the 
third  order  linear  dispersion  equation  on  a  periodic  domain, 
to  appear. 

Pa£e£s_b^  Dr_^  S^c£t^  ]![ • —Han^eti 

Frequency  -  Proportional  Damping  Mechanisms  in  Elastic  Beams, 
Doctoral  thesis.  University  of  Wisconsin-Madison ,  Dec.,  1988 

(with  G.  Weiss)  The  Operator  Carleson  Measure  Criterion  for  Ad¬ 
missibility  of  Control  Operators  for  Diagonal  Semigroups  on 
Systems  and  Control  Letters,  1991 

Pa£er.s_by.  P.ro.f^  ^ue_n^er_  Ite\age.ri,ng. 

A  Decomposition  Method  for  Integro-Partial  Differential  Equat¬ 
ions  and  Applications.  to  appear  in  Math.  Pures  et  Appl.,  1991 

Pa£e£s_bjr  Dr_^  ^egrg^e_Ve_is_Sj_ 
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Transfer  Functions  of  Regular  Systems:  Part  I,  submitted  to 
Trans,  of  the  AMS. 

Robustness  of  Feedback  Systems  with  Respect  to  Small  Time 
Delays,  Proc.  Conf.  on  Dec.  &  Contr.,  Honolulu,  Dec.,  1990 

Two  Conjectures  on  the  Admissibility  of  Control  Operators , 

Proc.  Vorau  Symposium  on  Ident.  &  Contr.  of  Distr.  Systems, 

July,  1990,  to  be  published  by  Birkhauser. 

Representations  of  Shift  Invariant  Operators  on  by  Transfer 

Functions  in  //“,•  An  Elementary  Proof,  a  Generalization  to 

and  a  Counterexample  for  t*.  Math,  for  Control,  Signals  and 
Systems,  1991 

5 .  Partially  Supported  Visitors  and  Speakers 

Speakers  partially  supported  by  the  grant  during  the  reporting  period 
with  the  general  area  of  their  presentation: 

Prof.  Katherine  Kime,  Case  -  Western  Reserve  University 

It 

Control  theory  of  the  Schrodinger  Equation 

Prof.  Joseph  Watkins,  University  oT  Southern  California 
Semigroups  and  Stochastic  F  -cesses 

Prof.  Luther  White,  University  of  Oklahoma 

Identification  of  Coefficients  in  Elliptic  Systems 


Prof.  B.  -  Y.  Zhang,  University  of  Cincinnati 

Control  Theory  of  the  Korteweg  -  de  Vries  Equation 

Prof.  Jack  Hale,  Georgia  Tech 

Periodic  solutions  of  infinite  dimensional  nonlinear  systems. 

Prof.  Walter  Rudin,  University  of  Wisconsin 

Some  aspects  of  the  theory  of  several  complex  variables. 
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Prof.  6.  Chen,  Texas  A.  4  M.  University 

Control  of  wave  equations  with  localized  controls. 

Prof.  W.  Hereman,  Colorado  School  of  .Mines 

Analytic  aethods  for  analysis  of  nonlinear  PDEs  with  solitary 
wave  Solutions 

Prof.  V.  Komornlk,  University  of  Bordeaux 

Minlaal  time  controllability  of  the  wave  equation  with  boun¬ 
dary  controls 

Prof.  L.  F.  Ho.  Wright  State  University 

Control  of  the  Euler  -  Bernoulli  beam  equation  with  localized 
controls 


6 .  Travel  to  Scientific  Meetings  Partially  Supported  by  AFOSR  89-0031 

Grant  funds  were  used  during  the  reporting  period  to  attend 
and  present  lectures  at  the  following  meetings.  In  each  case  a 
talk  was  presented  describing  AFOSR  -  supported  research. 

Conference  on  Decision  and  Control.  Tampa,  FL . ,  December,  1989 

SIAM  Meeting  on  Control  in  the  1990’s,  San  Francisco,  .May,  1989 

Three  Rivers  Applied  Mathematics  Symposium,  Pittsburgh,  PA, 

April  ,  1990 

SIA.M  Summer  Meeting,  Chicago,  Illinois,  June  1990 
Regional  AMS  Meeting,  Tampa,  Florida.  March,  1991 

(Travel  by  V.  Komornik)  Conference  on  Decision  and  Control,  Honolulu, 
December,  1990 

(Travel  by  G.  Weiss)  SIA.M  Meeting,  San  Francisco,  Feb.  1991 
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Soae  Reaarks  on  Experimental  Determination  of  Modal  Damping  Rates 

in  Elastic  Beams 
by 


David  L.  Russell 
Department  of  Mathematics 

Virginia  Polytechnic  Institute  and  State  University 

Blacksburg,  VA  24061 


1.  Background  Let  us  begin  with  a  protective  disclaimer  of  all  in¬ 
tention  to  treat  our  subject  in  any  degree  of  completeness,  particu¬ 
larly  in  regard  to  giving  an  account  of  historical  efforts  to  deter¬ 
mine  damping  rates  in  elastic  systems.  A  very  fine  account  of  the 
latter  is  provided  by  Bell  in  [  1  ] ,  particularly  in  regard  to  the 
nineteenth  century  studies  of  Kelvin.  Voigt  and  Weber.  An  overview 
of  recent  studies  in  this  direction  is  given  in  Bert's  article  [  2  ], 
We  should  also  duly  note  the  very  significant  study  [  12  ]  of  Zener 
in  which  he  presents  a  variety  of  theoretical  explanations  accounting 
for  internal  energy  losses  in  elastic  systems. 

The  projected  use  of  extended  elastic  structures  in  a  variety  of 
space  applications  has  led.  during  recent  decades,  to  a  significant 
revival  of  interest  in  the  study  of  internal  damping  mechanisms  in 
such  structures.  The  research  so  motivated  has  been  greatly  facili- 
ted  by  modern  methods  of  data  collection  and  analysis,  such  as  laser 
vibrometry  and  Fast  Fourier  Transform  (FFT)  techniques  allowing  far 
more  comprehensive  studies  of  damping  than  were  possible  in  earlier 
tines,  particularly  in  the  high  frequency  range.  Even  more  recently, 
as  illustrated  by  some  vibration  problems  encountered  with  the  Hubble 
space  telescope,  it  has  become  ever  more  evident  that  such  studies 
are  vital  if  future  space  platforms  are  to  function  satisfactorily. 

In  manned  space  stations  damping  is  likely  to  be  important  in  all 
frequency  ranges,  from  the  lower  frequencies  due  to  excitation  of 
natural  nodes  of  vibration  by  movement  of  masses  within  the  struc¬ 
ture,  docking  maneuvers,  etc.,  to  high  frequency  vibrations  excited 
by  rotating  electrical  machinery,  "popping"  of  joints  during  thermal 
expansion  or  contraction,  sharp  impacts,  and  many  other  causes. 
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Earthbound  laboratory  studies  of  Internal  damping  In  elastic 
structures  has  been  complicated  by  the  very  real  difficulty  of  eli¬ 
minating  energy  losses  to  the  external  environment,  which  are  often 
much  greater  than  the  internal  losses  which  are  the  subject  of  study. 
This  Is  particularly  true  in  the  high  frequency  range  where  vibration 
of  the  structure  excites  acoustic  waves  in  the  atmosphere  and  the 
support  platform  for  the  experiment.  Three  main  sources  of  external 
losses  may  be  cited:  i)  losses  due  to  attachment  of  sensing  devices 
such  as  accelerometers,  strain  gauges,  etc.,  li)  losses  due  to  exci¬ 
tation  of  sound  waves  in  the  atmosphere  and  atmospheric  viscosity, 
and  Hi)  losses  due  to  interaction  of  the  sample  with  Its  supporting 
structure,  such  as  a  clamp,  e.g..  The  first  of  these  can  be  largely 
eliminated  with  the  use  of  laser  vibrometer  measurements.  Ironically, 
the  use  of  that  technique  aggravates  the  problems  associated  with 
the  other  two  because  laser  vibrometry  is  a  sensitive  process:  it  is 
usually  difficult  to  obtain  a  satisfactory  reflected  signal  from  the 
sample  unless  the  sample  is  supported  in  such  a  way  as  to  minimize 
rotation.  This  problem  leads  to  the  necessity  of  using,  in  the  case 
of  samples  suspended  by  thin  nylon  filaments  ("pseudo-free"  configu- 
atlons),  fairly  large  samples  with  significant  rotational  inertia. 

In  the  case  of  lighter  samples  it  becomes  necessary  to  study  the  sam¬ 
ple  in  a  "cantilever"  configuration  using  a  rather  massive  clamp.  In 
either  case,  experience  indicates  significant  energy  losses  through 
the  stretched  filaments  or  the  supporting  clamp.  Laser  vibrometry  al¬ 
so  complicates  the  use  of  vacuum  chambers  to  reduce  atmospheric  eff¬ 
ects  since  a  laser  beam  often  cannot  be  transmitted  through  the  glass 
windows  of  such  a  chamber  without  significant  interference.  Thus  the 
vibrometer  needs  to  be  in  the  chamber  with  the  sample  with  attendant 
difficulties  in  aiming  the  laser  and  adjustment  of  the  associated  el¬ 
ectro-optic  devices. 

In  the  article  to  follow  we  present  the  results  of  the  writer’s 
experience  over  the  last  several  years  in  attempting  to  answer  two 
questions  about  internal  damping  in  elastic  beams.  These  questions 
are  the  following. 

1)  What  do  laboratory  experiments  indicate  concerning  the  qualita¬ 
tive  character  of  the  functional  relationship  between  modal  damping 
rates  arising  from  internal  energy  dissipation  and  the  corresponding 
modal  frequencies? 

ii)  What  precautions,  both  in  the  original  experiment  and  in  the 
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subsequent  Bathematlcal  analysis  of  the  experlnent  can  be  taken  to 
■iniaize  the  influence  of  external  energy  losses  so  that  the  observed 
daaping  rates  aay  reasonably  be  attributed  to  internal  dissipation? 

A  preliminary  exploration  of  both  of  these  questions  appears  in 
[  10  ] ;  it  amounts  to  an  Informal  account  of  our  understanding  of  the 
issues  and  the  implications  of  our  experimental  data  at  the  time  when 
that  article  was  written.  Developments  since  then  Indicate  that  our 
early  experiments  were  more  seriously  contaminated  by  the  effects  of 
external  losses  than  we  believed  at  that  time.  Also  in  the  interim 
a  better  understanding  has  been  gained  with  regard  to  the  mathemati¬ 
cal  modelling  of  both  Internal  and  external  energy  losses  in  elastic 
beams  and  of  mathematical  procedures  which  can  be  used,  in  conjunc¬ 
tion  with  corresponding  experimental  studies,  to  minimize  the  effects 
of  the  latter. 

In  §2  we  discuss,  from  a  mathematical  viewpoint,  what  losses 
may  be  expected  from  a  vibrating  beam  to  supporting  filaments  or  to 
the  surrounding  atmosphere.  Since  these  losses  appear  to  be  rela¬ 
tively  modest,  particularly  at  high  beam  vibration  frequencies,  we 
Indicate  at  the  end  of  that  article  the  nature  of  our  best  measure¬ 
ments  of  Internal  damping  in  certain  types  of  beams  suspended  by 
thin  nylon  filaments. 

In  the  article  [  10  ]  we  discussed  a  mathematical  approach  to 
determining  the  magnitude  of  energy  losses  from  a  beam  to  a  support¬ 
ing  clamp,  holding  the  beam  in  the  so-called  "cantilever"  configura¬ 
tion.  There  we  represented  the  clamping  device  as  a  much  thicker 
beam  to  which  the  thin  bean  under  study  is  coupled  in  a  particular 
way.  In  §3  we  present  a  revised  analysis  of  this  model  confirming 
the  results  obtained  in  a  more  acT  hoc  manner  in  that  article. 

Finally,  in  §4,  we  discuss  an  experimental  /  mathematical  proce¬ 
dure  based  on  observation  of  nodal  shapes  of  a  beam  subjected  to 
sinusoidal  forcing  which  shows  promise  of  being  able  to  identify  in¬ 
ternal  damping  rates  in  a  manner  minimally  affected  by  losses  to  a 
supporting  clamp.  We  exhibit  experimental  and  computational  data 
plots  to  indicate  the  promise  inherent  in  this  approach. 


-4- 


2 .  Losses  to  Supporting  Filaaents  and  the  Atmosphere 

As  we  have  Indicated  in  §1,  one  of  the  ways  in  which  one  may 
attempt  to  minimize  external  energy  losses  in  bean  damping  experi¬ 
ments  is  to  support  the  beam  by  long  thin  filaments,  commonly  nylon 
fishing  line,  as  shown  in  Figure  2.1  at  the  end  of  this  section.  In 
our  earlier  article  [  10  ]  we  provided  an  analysis  of  energy  losses 
from  the  beam  to  these  supporting  filaments.  Here  we  will  Indicate 
the  framework  used  in  that  analysis  and  restate  the  conclusions,  ref¬ 
erring  the  reader  to  the  original  paper  for  the  details  of  the  esti¬ 
mates  . 

Since  we  are  interested  in  energy  losses  from  the  elastic  bean 
to  the  supporting  filament,  we  will  use  the  energy  conservative  Euler 
-  Bernoulli  equation  to  model  the  beam.  Assuming  the  beam  to  have 
length  L,  uniform  linear  mass  density  n  and  bending  modulus  El,  all 
constant,  and  assuming  the  length  of  the  filament  supporting  the  beam 
to  be  t  and  the  linear  mass  density  to  be  r,  we  have  as  the  approxi¬ 
mate  energy  expression  for  the  system 


with  the  constraint  w(L,t)  »  v(0,t).  Here  w  is  the  bean  deflection 
as  a  function  of  its  longitudinal  parameter  x  and  v  is  the  filament 
deflection  as  a  function  of  its  longitudinal  parameter  s.  The  last 
term  in  the  first  line  cones  from  lifting  of  the  beam  as  it  swings 
on  the  supporting  filament.  From  this  energy  form  one  readily  derives 
the  equations  of  motion 


+  El 

^  -  0, 

0  X 

3t* 

ax* 

-  pL 

aif  =  0, 

0  <  s 

at* 

as* 

(2.01) 


(2.02) 


and  the  boundary  conditions 
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a*w 

ax* 


(o.t) 


0.  ^  (O.t) 

ax® 


0, 


a*w 

ax* 


(L.t)  =  0 


(L.t)  +  mL  ^  (O.t) 

ax®  as 


w  ( L  .  t ) 
t 


(2.03) 

(2.04) 


We  assume  that  waves  propagating  through  the  filament  are  not  reflec¬ 
ted  at  X  -  this  corresponds  to  the  boundary  condition 

^  (t.t)  +  ^(£.t)=0.  (2.05) 

at  as 

It  is  not  supposed  that  the  last  condition  is  completely  realistic. 

It  simply  replaces  more  complicated  assumptions  ensuring  that  very 
little,  if  any.  energy  entering  the  the  supporting  filament,  or 
string,  is  reflected  back  to  the  beam.  We  are  assuming  that  waves 
in  the  string  move  only  in  the  direction  from  s  =  0  to  s  =  t.  Then, 
in  fact,  equation  (2.05)  is  valid  with  t  replaced  by  s ,  0  <  s  <  £  ; 
in  particular  it  is  true  at  s  »  0.  Using  this  in  (2.04)  along  with 
the  constraint  w(L.t)  »  v(O.t)  we  arrive  at 

El  (L.t)  -  (rmL)^/*  ^  (L.t)  +  w(L.t)  .  (2.06) 

6x*  3t  ^ 

A  straightforward  spectral  analysis,  similar  to  that  appearing  in  the 
article  [  9  ].  shows  that  the  resulting  modal  damping  exponents  are 
asymptotically  constant  with  the  limiting  value 


It  is  eisy  to  see  that  the  attachment  of  n  parallel  filaments  corres¬ 
ponds  to  replacing  r  by  nr.  Thus,  as  we  vary  the  number  n  of  support¬ 
ing  filaments,  the  damping  rate  should  be  asymptotically  proportional 

1.  /  A 

to  n  .It  Is  clear  from  this  that  by  supporting  the  same  beam  in 
several  different  configurations  with  different  numbers  of  identical 
filaments  it  is.  at  least  in  principle,  possible  to  Identify  the 
energy  loss  rates  to  the  supporting  filaments  and  subtract  them  from 
the  overall  energy  loss  rates  to  provide  an  improved  estimate  of  the 
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Internal  dlssipatloa  in  the  bean. 

Let  us  now  turn  our  attention  to  the  question  of  losses  to  the 
surrounding  atmosphere  due  to  generation  of  acoustic  waves  from  the 
lateral  surfaces  of  the  beam.  It  must  be  recognized  that  there  may 
be  other  atmospheric  effects  than  these;  for  example,  viscosity  ef¬ 
fects.  particularly  due  to  airflow  around  the  edges  and  free  ends  of 
the  vibrating  beam.  These  nay  be  expected  to  have  very  complicated 
effects  and  we  make  no  attempt  to  account  for  these  in  the  present 
discussion.  In  order  to  keep  our  treatment  uncomplicated,  and  yet 
fairly  realistic,  let  us  consider  the  configuration  shown  in  Figure 
2.2;  the  ends  x  *  0  and  x  *  L  of  the  bean  are  identified,  so  that 
the  resulting  structure  is  L-periodic  in  the  x  direction.  For  high 
frequency  beam  vibrations,  which  are  our  main  interest,  acoustic 
waves  are  generated  predominantly  in  the  direction  orthogonal  to 
the  lateral  faces  of  the  beam.  Hence  we  study  the  linearized  acous¬ 
tic  equations  in  the  strip  S:  -»  <  y  <  «>  ,  0  <  x  <  L,  with  x  =  0 
again  Identified  with  x  >  L.  A  moment  of  consideration  leads  one  to 
conclude  that  atmospheric  density  variations  induced  by  motion  of 
the  beam  will  be  odd  functions  of  y;  hence  we  may  confine  our  study 
to  the  half  strip  corresponding  to  y  >  0.  We  begin  with  the  Euler 
equations  for  compressible  flow  in  two  dimensions,  as  given  in  [  6  ] 


ff . ) , 

for 

example : 

p 

au 

+  P 

u  ^  + 

p 

V 

1“  + 

P'  ^  =  0  , 

(2.07) 

at 

ax 

3y 

ax 

p 

H 

+  P 

u  + 

p 

V 

P-  ^  =  0  . 

(2.08) 

at 

ax 

ay 

3y 

i£ 

+  u 

a£ 

X  V  9l£ 

+ 

p 

f  ^ 

^^1=0. 

(2.09) 

at 

ax 

ay 

ax 

ax  ^ 

where  p  is  the  atmospheric  density,  u.v  are,  respectively,  the  fluid 
velocity  components  in  the  x,y  directions,  and  p  »  p(p)  describes 
the  constitutive  relation  between  pressure  and  density.  To  obtain 
linearized  equations  for  small  amplitude  waves  we  let  be  the  mean 

atmospheric  density  and  we  represent  p  in  the  form  p  =  p ^  l  6p, 

u  =  t  p .  V  *  e  1/ .  Setting  ft  =  p'(P|j)  and  collecting  terms  of  first 

order  in  e  in  (2.07)  -  (2.09),  finally  multiplying  the  last  equation 
by  ft/p^,  we  arrive  at  the  linearized  equations 
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p  iii  +  R  iij2  =  0  , 

°  at  ax 

(2.10) 

P„  ^  R  ^  =  0  . 

°  at  ay 

(2.11) 

(ff/p„)  ^  +  Jt  f  ^  ^  I  =  0  . 

“at  '■ax  ax 

(2.12) 

This  is  a  linear  symmetric  hyperbolic  system  with  conserved  energy 
integral 


g  - 


^0^ 


]  dx  dy  . 


(2.13) 


defined  (minimally)  on  any  bounded  region  R  contained  in  the  half 
strip  described  above. 


Almost  any  linear  operation  on  the  variables  n,  v,  6p  will  pro¬ 
duce  a  solution  z  of  the  wave  equation 


at* 


f  3^  +  liz  ] 
^  ax*  ay*  J 


(2.14) 


the  trick  is  to  find  the  linear  operation  for  , which  the  energy  asso¬ 
ciated  with  (2.14),  i.e., 


E  - 


[  t;  ]’]}  - 


(2.15) 


agrees  with  the  physical  energy  (2.13).  Assuming  the  mean  value  of 
each  of  p,  v.  6p  to  be  zero,  we  may  define  the  inverse  Laplace  opera¬ 
tor  A~^  on  each  of  these  quantities  or  their  derivatives.  Applying 
a*/atax  to  (2.10)  and  3*/atay  to  (2.11)  and  adding,  we  have 


p  -div(;*,  1/) 

”  at* 


Jf 


— A«5p 

at 


0 
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Applying  we  have,  with 


z  =  A  ^div(p,v), 


(2.16) 


a^z 

at* 


R  ^ 

at 


•  ill 
°  at* 


-  ^0^ 


div(M.v)  =  -  Jt  AZ  ]  = 


and  we  conclude  that  (2.13)  is  satisfied  with  z  given  by  (2.16). 
Applying  a/ax  to  (2.10)  and  a/ay  to  (2.11),  adding  and  then  applying 
A~^ ,  we  have 


az  R 

at  ^0 


(2.17) 


Using  this  equation  to  solve  for  5p  in  (2.10),  (2.11)  we  have 


f  M  -  -  1 

=  0 

- 

'  u  - 

=  0 

at 

^  ax  J 

at 

ay  ^ 

From  this  it  can  be  shown,  auch  as  in  [  8  ] ,  that 

az 


ax 


3y 


Then  (2.13)  becomes 


*  ■  »  n  fy  If  ^ 


and  we  see  that  z.  as  given  by  (2.16),  yields  equivalent  energy  int¬ 
egrals;  indeed,  replacing  z  by  the  energy  integrals  agree 

Accordingly  the  total  energy  expression  for  the  coupled  beam  - 
atmosphere  system  becomes 
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Dif f erentiating  with  respect  to  t,  assuming  the  overall  system  to  be 
energy  conservative,  and  recognizing  that  the  force  due  to  the  pres¬ 
sure  variation  across  the  beam  leads  to  the  equation 


m 


a^w 

at® 


El 


a*w 

ax* 


2  e  6  p 


0  . 


the  usual  calculations  lead  to  the  interface  condition 


pL 

=  [c  ^(x,t)  6p(x,0.t)  +  p  ^(x,0,t)  -^(x.O.t)] 

Jq  L  at  °  ay  at  J 

pL 

\t  ^(x,t)  -  ff  ^{x  3,t)l  ap(x.o.t)  . 

Jn  ‘■at  ay  j 


dx 


from  which  we  conclude  that  the  appropriate  boundary  condition  for  z 
is 


^(x.O.t)  =  *  ^(x.t)  .  (2.18) 

ay  at 

Let  us,  for  convenience,  normalize  L  to  2r.  We  introduce  the 
Fourier  expansions 


w(x, t) 


=  i  w.  (t)e^‘^*, 
k  =  -oo 


z(x,y,t) 


=  i  zj^{y,t)e^^'', 

k=-oo 


ap(x,y,t)  =  2  (y . t ) e^^^  . 

k=-«  ^ 

Then,  with  «  =  El/m,  the  coefficients  ( t )  satisfy  the  ordinary 
differential  equations 


d^ 

dt* 


-k  *  a  Wj^  (  t  ) 


2e  <5pj^(0,t),  -«  <  k  <  o® 


(2.19) 


while  the  coefficients  Zj^(y,t)  obey  the  partial  differential  equa¬ 
tions 
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3*2j 

3t* 


-  R 


3*Zj 

3y® 


R  k“  Zv  =  0  .  <  k  <  oo  , 


(2.20) 


and  the  boundary  conditions  at  y  =  0  (cf.(2.18)) 


3zi 


-(o.t)  =  i  — Ht)  ,  -»  <  k  < 
3y  dt 


Using  (2.17)  we  may  replace  (2.19)  by 


(2.21) 


d*w 


dt‘ 


^+k*awv(t)  - 


2tp^  azjc 

~ir  Ti- 


(O.t)  =0,  -~<k< 


(2.22) 


In  order  to  be  able  to  use  (2.18)  to  replace  ^^k  in  (2.22)  by  *^”k  , 

Ft“ 

we  need  to  establish  a  connection  between  the  t  and  y  derivatives 
of  Zj^  at  y  =  0.  To  this  end  we  need  a  more  detailed  description  of 

the  forced  oscillations  in  the  atmosphere.  To  this  end  we  attempt 
a  solution  of  (2.20)  -  (2.22)  in  the  form 


zj^(y.t)  =  fk  ®  "k  ® 


(2.23) 


Then  from  (2.20)  we  see  quite  directly  that 

X*  =  K  -  k*  , 

while  (2.21)  and  (2.22)  combine  to  give 

+  ak  * 

^K-^Xj^ 


(2.24) 


Setting  the  determinant  of  the  matrix  equal  to  zero  and  substituting 
('5  24),  we  have  the  polynomial  equation  in 
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K  -  2cPgK~‘  p*  +  (ak*  -  k*)  +  2e*P|jK'®k*  =  0  .  (2.25) 

Let 

Pj,  =  k®  .  (2.26) 

<w 

Substituting  (2.26)  into  (2.25)  and  dividing  by  Mj^k®,  we  have,  as 
k  ->  <»  , 

K  mJ  -  2c*PjjK"^k"*Pjj  +  a  -  0(k”*)  =  0  . 

which  gives 

Mk  *  c*PgK~*k~*  ±  i  («K"M‘''*  +  <5(k“*) 

and  thus 

Mjj  »  ±  i  («K"‘)‘/*k*  +  «*PoK"*  +  0(k"*)  .  (2.27) 

Substituting  (2.27)  into  (2.24)  we  find  that 

Xk/Mj^  =  ±  +  0(k"*).  k  0®  .  (2.28) 

The  -  sign  is  used  since  it  corresponds  to  waves  moving  away  from 
the  elastic  beam  in  the  positive  y  direction.  Then  an  easy  compu¬ 
tation  using  the  assumed  form  (2.23)  and  (2.28)  shows  that 

®®k  (0,t)  =  f-K*/®  -  (3(k”®)l  ®®k  (0,t)  =  f-tK"^/®  +  (3(k'®)]  f![k  (t) 

TT*  •*  W~  ■'  cTt^ 

and  thus  (2.22)  becomes 

^  r2e*p  K”®'^®  +  0(k“®)l  ‘^^k  +  k*  a  w^  ( t )  .  -«>  <  k  <  <»  . 

dt®  ^  J  HT“ 

From  these  equations  we  conclude  that  rate  of  energy  dissipation  in 
the  k-th  vibrational  mode  is  asymptotically  uniform  as  k  ->  »  and 
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dlrectly  proportional  to  the  sean  atmospheric  density  p^. 

In  Figure  2.3  we  show  superimposed  plots  of  successive  power 
spectra  for  an  aluminum  bar  supported  by  two  nylon  filaments  at  nor¬ 
mal  atmospheric  pressure;  the  delay  from  the  beginning  of  the  first 
sampled  time  Interval  to  the  beginning  of  the  second  is  2.8  seconds. 
The  lighter  curve  is  the  power  spectrum  corresponding  to  the  first 
interval;  its  peaks  are  highlighted  with  a  dark  dot.  The  darker 
curve  corresponds  to  the  second  power  spectrum.  (See  [  9  ]  for  a 
more  detailed  account  of  the  manner  In  which  these  successive  power 
spectra  are  obtained.)  Figure  2.4  plots  the  differences  of  the  peak 
values,  corresponding  to  the  damping  rate  in  db/sec  for  the  natural 
modes  of  vibration  with  the  Indicated  frequencies.  These  plots  very 
strongly  indicate  a  rate  of  damping  at  least  linearly  proportional 
to  frequency.  The  anomalously  high  damping  rates  at  the  lower  end 
of  the  spectrum  are  due  to  losses  to  the  resonant  modes  of  the  sup¬ 
porting  filaments.  Our  mathematical  analysis  Indicates  that  we  need 
not  anticipate  that  the  linear  character  of  the  damping  versus  freq¬ 
uency  relationship  in  the  high  frequency  range  is  seriously  corrupted 
by  losses  to  the  supporting  filaments  or  the  atmosphere. 


Damping  Exponent 


3 .  External  Energy  Losses  for  Clamped  Beams 

As  in  our  earlier  article  [  10  ].  we  consider  a  thin  beam  whose 
lateral  displacement  is  described  by  w(x.t),  0  <  x  <  L.  This  beam  is 
assumed  clamped  at  x  >  0  by  a  device  which  we  model  mathematically  as 
a  thick  beam  with  displacement  W(x.t).  We  consider  "monochromatic", 
i.e.,  single  frequency,  solutions  of  the  combined  system  under  the 
assumption  that  "outgoing"  waves,  moving  toward  -<»  in  the  clamping 
beam  are  the  only  motion  present  in  that  structure.  For  definiteness 
we  assume  the  x  =  L  endpoint  of  the  thin  beam  to  be  free  but  other 
energy-conservative  boundary  assumptions  would  serve  equally  well. 

Our  analysis  here  will  differ  somewhat  from  that  given  in  [  10  ]  , 
where  we  assumed  the  thin  beam  to  have  infinite  length.  We  refer 
the  reader  to  that  article  for  explanatory  figures. 

We  treat  the  thin  beam  in  the  Euler  -  Bernoulli  framework  with 
equation  of  motion 

m  -  el  =  0  ,  (3.01) 

at*  ax* 

a  being  the  mass  density  per  unit  length  and  ei  the  bending  modulus. 
At  X  »  L  we  impose  the  familiar  free  endpoint  boundary  conditions. 

The  thick  beam  to  which  the  thin  beam  is  attached  is  modelled  via  a 
similar  equation  with  corresponding  parameters  M  and  EI; 

M  Sll  ^  gj  a^  ^  Q  (3.02) 

at*  ax* 

The  clamping  assumption  yields  the  kinematic,  or  essential,  boundary 
conditions 


w(O.t)  ■  W(0,t) . 

The  relevant  energy  expression  is 


3W(0,t)  =  ^(O.t) 
ax  ax 


(3.03) 


(3.04) 


and  it  is  easy  to  see  that  the  conditions  which  should  be  satisified 
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for 

energy 

to  be 

conserved  at  the 

Joint  X  a  0  between  the 

thin 

and 

thick  beams 

is 

El 

^(0,t) 

■ 

ei 

^(0,t)  , 

(3.05) 

8x® 

El 

^(O.t) 

B 

el 

^(0,t)  . 

(3.06) 

3x* 

8x* 

Let  us 

begin 

the 

analysis 

by 

studying  solutions  of  (3 

.02) 

with 

the 

special 

form 

W(x.t) 

e^*  . 

(3.07) 

We  further  restrict  these  solutions  by  insisting  that  they  be  "out¬ 
going  solutions",  i.e.,  that  they  should  represent  waves  moving  away 
from  X  -  0  toward  -<».  To  facilitate  analysis  of  what  that  means,  let 
us  set 


X  3  Clfl 


(3.08) 


where 


a  ■  1  (3.09) 

comes  from  a  set  of  complex  numbers  yet  to  be  determined.  Substitu¬ 
ting  (3.07),  (3.08)  into  (3.02)  and  dividing  by  we  arrive  at  the 

equation 


M  a*  +  El  M*  -  0 


so  that 


M  =  ±  v' ^ / E 1  ai  «  ±  R  ai  .  (3.10) 

Using  (3.10)  and  specializing  to  the  +  sign  for  the  present,  we  have 


Ri  [(<*;-«»)  t 


«oX] 


[-2«  «  t-a  x] 

°  ^  ^  (3.11) 


W( X , t )  =  e 


R1 ( a*t  +  «x ) 
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The  condition  for  waves  to  be  noving  fron  right  to  left  is  that 

and  should  have  the  sane  sign.  Since  this  notion  is  supposed 

to  result  fron  a  notion  in  the  thin  bean  decaying  in  time,  both  a  a 

0  1 

and  in  the  real  exponent  should  be  positive.  Thus  we  are  looking 

for  positive  nunbers  and  a  with  <x  >  a  . 

0  10  1 

The  corresponding  displacement  w(x,t)  of  the  thin  beam  will  be 
assumed  to  have  the  form  (cf.  (3.01)) 

w(x,t)  “  e^^<p(x)  =  e^*  ^  <p(x)  ;  (3.12) 

substitution  into  (3.01)  provides  the  differential  equation  satisfied 
by  ? (x) : 


ei  ^  =  ma*9  , 
dx* 

so  that,  with  r  *  (n/el)^'^*, 

«(x)  «  c^cos  rax  +  c^sin  rax  +  +  c  ,  x  c  (0,L). 

‘  ®  *  (3.13) 

Fron  the  formula  (3.07)  for  W  and  the  identities  (3.08),  (3.10), 

we  have 


3^W 

3t3x 

(x .  t ) 

s  X/i 

e’^^e^''  »  -Ra®  e^^e^* 

3W 

3t 

(x ,  t ) 

*  X  e 

=  Ra*i  e’^^e^** 

EI 

3*W 

3x* 

(X,  t) 

=•  -  EI  R  a*  e^^e^’'  , 

EI 

3^W 

3x* 

(x,t) 

-  -EI  R®a®i  e^^e^*  . 

Then  fron  (3.05),  (3.06)  and  the  requirement  that  the  boundary  condi¬ 

tions  on  w  at  X  >>  0  should  reduce  to  the  free  endpoint  conditions 
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as  M  and  El  tend  to  0.  we  obtain 


el  ^(O.t)  -  El  ^(O.t)  •  -  El  R*«  ^(O.t)  .  (3.14) 

8x*  ax®  at 


el  ^(O.t)  «  El  ^(O.t)  -  El  ^  (O.t)  .  (3.15) 

ax*  ax*  atax 

When  the  boundary  conditions  (3.14).  (3.15)  are  applied  to  coaplex 
solutions  of  (3.01)  along  with  the  conservative  free  endpoint  condi¬ 
tions 


one  finds  that  the  tine  rate  of  change  of 


E 


(3.16) 


Is 


-Re  f  el  ^(O.t)  ^  (O.t)  -  el  ^(O.t)  i^i(O.t)  ]  - 

L  ax*  atax  ax*  at  ■> 

-Re  f  el  ^(O.t)  ^  (O.t)  -  el  3l^(0.t)  ^(O.t)  1  * 

L  ax*  atax  ax*  at  J 

-El  Re(a"‘)  ll-i  (O.t)!  -  El  R*  Re(a)  |ill.(0,t)|  <  0  (3.17) 

•atax  I  lat  ' 


with  our  assuaptlon  that  Re(«)  -  >  0.  The  remainder  of  the  analy¬ 
sis  consists  of  an  examination  of  the  system  (3.01),  (3.14)  -  (3.16) 

with  a  view  to  obtaining  the  high  frequency  nodal  damping  rates. 

From  the  boundary  conditions  (3.14),  (3.15),  (3.16)  and  from 

(3.12) .  setting 


r  =•  El/ei ,  4  =«  (El  R*)  /ei  , 
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we  obtain  for  (p  the  boundary  conditions 

<p'  '  (L)  «  0.  <p'  ’  '  (L)  =  0.  (3.18) 

<p''(0)  =  jai  <p'(0).  <p''’(0)  =  -«a®i  <?  ( 0 )  .  (3.19) 

Using  the  so-called  wave  propagation  method  [  4  ],  which  simplifies 
eigenvalue  approximations  by  neglecting  the  exponentially  small  terms 

g-raL  ai.i3j[ng  from  substitution  of  (3.12)  into  (3.18),  (3.19)  we  ob¬ 

tain  the  vector  equation 

•  -cos  r«L  -sin  raL 
sin  raL  -cos  raL 
-r  -yi 


Setting  the  determinant  of  the  indicated  matrix  equal  to  zero  we  have 
sin  raL  ^(ayr*  +  2r8)lj  +  cos  raL  |^2r*  +  2y6  +  {2yT^  -  28r)ij  »  0, 
or 

tan  r«L  =  2Ar..  -  27..^  t-(2r*  ZySU 

2rr®  +  2r<5 

Since  tan  raL  behaves  like  l/(  (k  +  l/2)TT-raL)  in  the  neighborhood  of 
(k  +  l/2)Tr  for  each  Integer  k  we  have,  in  first  approximation 


1  0 

1  0 

0  7l  +  r 
-r®+8 i 


(3.20) 


raL 


(k+l/2)r 


_ 2Tr^  +  2t6 _ 

2«r  -  2yr^  +  (2r*  +  27<5)i 


-  (k  +  l/2)r  -  (2rr®  -  2r«)  ~  ^7^^)  -  .(^r*  +  27<5)i 

(2«r  -  2rr*)®-^  (2r*  +  27«) 

and  we  conclude  that  large  values  of  a  do,  indeed,  have  the  form 

(3.09)  with  a^  and  a^  positive;  a^  grows  like  (k+l/2)T:/rL  +  c  for 

some  real  constant  c  while  a^  tends  to  a  fixed  positive  value.  Ac¬ 
cordingly  , 
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X  »  Ra®i  *  ^[(*0  ~  “  ”^^*0*1  '*’  ~ 

has  negative  real  part  «  -2Ra^«^  asyaptotlcally  proportional  to  the 

square  root  of  the  frequency  a  R(«q  ~  «* )  as  the  latter  tends  to 

+  «>.  The  other  values  of  X  being  complex  conjugates  of  these,  we 
conclude  that,  granted  the  validity  of  this  type  of  clamping  model, 
the  damping  rate  should  be  asymptotically  proportional  to  the  square 
root  of  the  frequency.  This  agrees  with  the  result  obtained  in  [  10  ] 
via  slightly  different  reasoning. 


There  may  be  some  disagreement  that  the  thin  beam/thick  beam 
representation  of  a  clamping  device  is  an  accurate  one.  The  assump¬ 
tion  of  infinite  length  for  the  thick  beam  is  certainly  suspect  as  is 
the  assumption  that  no  energy  is  lost  in  the  clamped  joint  Itself. 

A  somewhat  different  model  is  obtained  if  we  abandon  the  thick  beam 
representation  of  the  clamp  and  suppose  the  bean  attached  at  x  »  0 
to  a  clamping  device  which  may  undergo  some  rotational  deformation  in 
the  process  of  bean  oscillation  and  nay  dissipate  energy  in  that  pro¬ 
cess.  For  this  model  we  suppose  an  overall  energy  functional 


The  last  tern  corresponds  to  rotational  elasticity  of  the  clamp;  if 
the  clamp  were  assumed  perfectly  rigid  one  would  impose  a  boundary 
condition  3w/3x(0,t)  =  0.  For  smooth  solutions  we  nay  then  compute, 
taking  the  boundary  conditions  at  x  =  L  to  be  (3.16)  and  employing 
the  usual  integrations  by  parts  in  the  second  term  of  the  Integral, 


^  -  El  ^(0,t)-^(0,t)  +  EI^(0,t)^(0,t)  +  R^(0,t)  ^(0,t). 

dt  3x*  3t3x  3x®  3t  3x  3t3x 


Assuming 


w(0 , t)  mo. 


(3.22) 


the  first  and  third  terms  reduce  to 
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^^(O.t)  Tr  ^(O.t)  -  El  ^(O.t)]  . 

3t3x  L  Zx  ax*  J 

which  becomes  the  boundary  dissipation  tern 

2  2 

-a  r-^-(O.t)]  <  0  (3.23) 

•■atax  J 

if  we  impose  at  x  >  0  a  second  boundary  condition 

R  ^(O.t)  -  El  ^(O.t)  +  a  -^{O.t)  =  0  .  (3.24) 

ax  ax*  atax 

wherein  the  last  term  corresponds  to  an  assumed  frictional  resist¬ 
ance  to  rotation  in  the  supporting  clamp. 

Analysis  of  the  system  consisting  of  the  equation  (3.01)  along 
with  the  boundary  conditions  (3.16),  (3.22)  and  (3.24)  by  the  wave 

propagation  method  referred  to  earlier  shows  that  in  this  case  the 
rate  of  energy  decay  in  each  mode  is  proportional  to  the  modal  fre¬ 
quency.  If  this  model,  or  even  the  model  discussed  above,  for  energy 
loss  to  a  supporting  clamp  is  a  valid  one  it  is  clear  that  the  deter¬ 
mination  of  Internal  damping  rates  in  clamped  beams  from  Just  the 
"raw"  spectral  data  is  a  questionable  undertaking.  In  the  final  sec¬ 
tion  of  this  article  we  indicate  mathematical  procedures  whereby  the 
clamping  effects  can  be  very  nearly  eliminated  through  observation  of 
modal  shapes  in  situations  where  the  supporting  clamp,  and  thus  the 
beam  with  it,  can  be  forced  to  oscillate  at  the  corresponding  modal 
frequency . 

In  view  of  the  above  results  it  is  significant  that  the  writer 
has  never  been  able  to  obtain  satisfactory  damping  data  from  clamped 
metallic  beams  in  the  laboratory,  notwithstanding  any  earlier  conclu¬ 
sions  in  [  10  ] .  Ne  have  obtained  fairly  good  data  for  composite 
beams  under  these  circumstances  because  they  have  such  high  internal 
damping  rates  that  those  rates  are  not  unduly  masked  by  losses  to  the 
clamping  device. 
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4 .  Characteristic  Features  of  Elastic  Beams  In  Forced  Oscillatory 
Motion  and  their  Relation  to  Internal  Damping  Rates 

We  provide  here  a  mathematical  analysis  showing  that  an  elastic 
beam  excited  into  vibration  by  sinusoidal  motion  of  its  clamping 
mechanism,  exhibits  in  the  resulting  beam  vibration  certain  features 
which  are  determined  by  the  Internal  damping  rate  and  are  largely 
independent  of  the  characteristics  of  the  clamp.  In  Figure  4.1  at 
the  end  of  this  section  we  indicate,  schematically,  the  experimental 
set-up  and  the  heuristic  basis  for  the  proposed  method. 

The  forced  motion  of  the  clamped  end  of  the  beam  generates  waves 
propagating  "outward"  toward  the  free  end  of  the  beam,  where  they  are 
reflected  and  travel  back  toward  the  clamped  end.  In  making  this  dou¬ 
ble  transit  of  the  beam  the  waves  are,  to  a  greater  or  lesser  extent, 
depending  on  the  degree  of  internal  damping,  attenuated;  the  work 
done  by  the  forcing  function  at  the  clamped  end  is,  in  the  periodic 
steady  state  situation,  just  enough  to  make  up  for  the  energy  lost 
in  the  propagation  process  together  with  losses  to  the  supporting 
structure  . 

In  the  loss-free  situation,  once  a  periodic  motion  has  been  in¬ 
itiated,  it  continues  without  any  motion  of  the  clamped  end.  The 
waves  returning  from  the  free  end  have  the  same  amplitude  as  the  out¬ 
going  waves  and,  at  the  natural  frequencies  of  the  beam,  the  outgoing 
and  incoming  waves  interfere  constructively  to  produce  a  standing 
wave  which  is  the  corresponding  natural  mode  of  vibration.  When  we 
have  internal  damping,  however,  the  returning  wave  has  smaller  ampli¬ 
tude  than  the  outgoing  wave  and  the  constructive  interference  is  not 
perfect.  In  particular,  as  we  will  see,  there  are  no  true  "nodes" 
and  the  steady  state  amplitude  is  not  the  absolute  value  of  a  sinus- 
soidal  function  of  x  (plus  the  small  correction  due  to  the  real  expo¬ 
nential  parts  of  the  eigenfunction  in  question)  but,  rather,  is  modi¬ 
fied  so  as  to  decrease  in  amplitude  as  x  varies  from  the  clamped  end, 
X  =  0,  to  the  free  end,  x  *  L.  We  will  see  that  this  spatial  varia¬ 
tion  in  amplitude  provides  a  measure  of  the  internal  damping  in  the 
beau  which  is  independent  of  energy  losses  in  the  clamping/forcing 
structure  holding  the  beam  at  x  =  0. 

In  our  present  mathematical  discussion  we  will  assume  that  the 
Internal  damping  operator  commutes  with  the  elasticity  operator.  This 
is  not  the  case  for  many  recent  beam  dissipation  models  (cf.[  9,10  ]) 
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but  it  is  true  for  viscous  damping,  Kelvin-Voigt  damping  (cf.  [  10  ]) 
and  square-root  damping  [  3,5,7  ]:  it  will  serve  for  the  present  ana¬ 
lysis.  Accordingly  we  consider  a  uniform  elastic  beam  modelled  by  an 
an  Euler  -  Bernoulli  type  equation 

P  +  G  ^  +  El  =  0,  t  e  (-~,«)  ,  X  6  [0,L]  ,  (4.01) 

Bt*  at  ax* 

where  G  is  a  non-negative,  in  general  unbounded,  self-adjoint  opera¬ 
tor  commuting  with  the  operator  a*/ax*.  The  beam  is  assumed  free  at 
the  end  x  *  L,  corresponding  to  boundary  conditions  (3.16).  At  x  =  0 
we  suppose  that  there  is  forced  motion: 

w(0,t)  =  f(t)  .  (4.02) 

This  motion  is  produced  by  a  corresponding  motion  of  a  clamping  de¬ 
vice.  We  model  the  clamp  as  in  the  last  paragraphs  of  the  preceding 
section,  supposing  the  overall  energy  functional  to  be  (3.21),  The 
rate  of  change  of  energy  is  now  seen  to  be 


dS  =  _  r3W(.,t),  6  ^(•,t)l  -  El  ^(O.t)  -^(O.t) 

dt  '•at  at  "^L^fo  L]  atax 


+  El  ^^(O.t)  ^(0,t)  +  R  ^(O.t)  -^^(0,t)  .  (4.03) 

3x^  at  ax  atax 

The  first  term  on  the  right  hand  side  is  non-positive  and  represents 
internal  energy  dissipation.  From  (1.03)  the  third  term  is 


El  ^(0,t)  f(t)  ,  (4.04) 

ax® 

representing  the  work  done  on  the  beam  by  the  external  forcing.  With 
the  boundary  conditions  (3  t)  and  (4.02)  applying  at  x  =  0  the  total 
energy  dissipation  form  is  ow  seen  to  be 


L®[0,L] 


(4.05) 


which  is  balanced  by  (4.04)  over  each  period  of  the  steady  state  for- 
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ced  oscillation. 

Now  we  suppose  that  the  function  f(t)  in  (4.02)  takes  the  form 

f(t)  -  e^“^  (4.06) 

and  that  the  corresponding  response  takes  the  form 

w(x,t)  =■  w(x,u)  :  (4.07) 

we  will  separate  out  the  real  component  of  the  motion  at  the  end  of 

our  calculations.  The  boundary  conditions  for  w  now  take  the  form 
(cf.(4.02),  (3.16),  (3.24)) 

w(0.«)  ■  1  .  (R  -  al<y)  ^(O.m)  -  El  ^^(0,«)  ■  0  .  (4.08) 

3x  3x* 

A  ^ 

^(L,«)  ■  ^(L,«)  «  0  ,  (4.09) 

3x*  3x® 

From  our  commutativity  and  non-negativity  assumptions  on  the  operator 
G,  we  may  assume  that  there  is  a  non-negative  function  7(«),  the  form 

of  which  is  to  be  Identified,  such  that  w(x,<i;)  satisfies  the  differ¬ 
ential  equation 

El  ^(x,<a)  +  (iuy(u)  -  p«*)  w(x,«)  =  0  .  (4.10) 

dx* 

Suppressing  the  role  of  «  as  an  argument  and  setting 

,  ^  4i(TT/2-E) 

r  =  p/EI,  «g  =  y/El,  w*(r-ig)  =  v*e  ,  (4.11) 

where 


41  (7r/2-£) 


=  (r-ig)  / 


/./;t 


the  fourth  roots  of  «*(r-ig)  take  the  form 
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1  ( kTr/2-c  ) 

Pj^»i/e  »i/ej^,k  =  l,2.3,4.  (4.12) 

The  Indexing  is  such  that  lies  in  the  k-th  quadrant  of  the  complex 
plane  in  each  case,  assuming  g  sufficiently  small  so  that  0  <  e  <  . 

Accordingly,  for  some  coefficients  Cj^,  k  =  1,2, 3, 4,  we  have 


fi  X  p  X  p,X  ft  X 

w(x,m)  =c^e^  +ce*  +ce^  +  c  e  (4.13) 

1  2  3  V 

Applying  the  boundary  conditions  (4.08),  (4.09)  we  have  the  system  of 

equations 


C,  ^  +  Cg  +  C^  =  1  , 


(4.14) 


(4.15) 


2  ^1^  a  a  a  ^4^ 

u  e  c  +  p  e  c  +  u  e  c  +  u  e  c  =  0  , 

'^1  ^2  a  ^3  3  ^4  4  ' 


(4.16) 


a  a  ^2^  a  ^a^  a  ^4^ 

p®e  c^  +  fi\e  *  +  p|e  *  c^  +  p®e  *  c^  =  0  , 


(4.17) 


where,  in  agreement  with  (4.08), 


Rj^  =  (R  +  ol«)  -  El  pj  ,  k  =  1,2, 3, 4  . 

The  solution  of  this  four  dimensional  system  is  rather  compli¬ 
cated;  it  will  serve  our  purposes  here  to  develop  an  approximate  sol¬ 
ution,  valid  for  large  v  and  small  z,  again  using  the  wave  propaga¬ 
tion  method.  Under  these  circumstances  p^  has  large  positive  real 

part  and  p  has  large  negative  real  part.  Assuming  c  and  c  can  be 

shown  to  have  relatively  moderate  values,  we  omit  the  second  term 
from  equations  (4.16)  and  (4.17).  .Multiplying  the  modified  equation 

(4.16)  by  ft ^  and  subtracting  the  modified  equation  (4.17)  we  have 


(4.18) 


-26- 


Fro«  (4.18)  we  Infer  that  must  be  very  small.  Omitting  the  fourth 

term  from  equations  (4.14),  (4.15),  we  obtain  two  modified  equations. 
Then,  multiplying  the  modified  equation  (4.16)  by  and  subtracting 

the  modified  equation  (4.17)  we  obtain  a  third  equation.  These  toget¬ 
her  form  a  three  dimensional  system 

+  Cg  =  1  ,  (4.19) 

^3*^3  ”  °  ’  (4.20) 

e^^  (m^M*  -  M®)c^  +  e^®  "  °  •  (4.21) 

Solving  for  in  (4.19)  and  substituting  in  (4.20)  and  (4.21)  we 
obtain  the  two  dimensional  system  consisting  of  (4.21)  and 

(R2-Ri)c^  +  (^2'‘^3^'^3  '  *^2  •  (4.22) 

Applying  Cramer's  rule  to  (4 . 21 ) , ( 4 . 22 )  we  have 

M  b 

=  R^e  ®  (m^M®  -  M®)/C(M)  .  (4.23) 

M  b 

Cg  =  -R^e  ‘  ~  ,  (4.24) 

where  U(m)  is  the  determinant 

2J(m)  =  "  ^3^  "  .(4.25) 

From  (4.23),  (4.24),  noting  that  /*  »  -  p  ,  we  see  that  there  is  a 

3  i 

complex  A  =  A(R,o,m)  such  that  (cf.(4.12)) 

M,b  M,b 

Cl  =  A  (e^  -  eg)  e  ®  ,  Cg  =  A  (ei  -  e^)  e  ^  (4.26) 

From  (4.18)  we  also  have,  since  m*  = 
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c 


4 


( M  I* 


e  * 


(4.27) 


IQ  this  approximation  e  =  e  is  to  be  regarded  as  replaceable 

P.L 

by  zero,  unless  multiplied  by  a  term  comparable  to  e  .  Thus  we 
have ,  using  (4.19) , 


1  -  c,  -  c. 


1  -  A 


p  L 

e^  (e  + 


-M,  L 

e  "  )  - 


M,  L 

e^  (e 


-  e 


>] 


(4.28) 


Thus  the  forced  beam  oscillation  has  the  approximate  form 


- 1  »■  <  r  (x-L)  -n  IX.-L,)  . 

w(x,t)  =  Ae  ^  e^  ^r(l+l)e  ^  +(l-l)e  ^  +(l-i)e^  e 


-M, (x-L) 


ic  m^(x-L) 


+  e 


^®[a  ‘  -  [e^(e  ‘  +  e  *  )  -  e^(e  ^  -  e  ^  )]]e  *  j  ,(4.29) 


Since  A  depends  on  a.  the  damping  parameter  of  the  clamp,  the  last 

-M.x 

term  also  depends  on  a.  But,  since  e  is  a  factor  in  this  term 

and  since  has  large  negative  real  part  for  large  values  of  u, 

this  term  will  not  be  important,  again  for  large  u,  at  any  substan¬ 
tial  distance  from  the  clamp.  The  third  term  will  not  be  important 
at  any  substantial  distance  from  the  tip  of  the  beam.  Assuming  we 
make  measurements  in  a  "central  region"  of  the  beam,  we  can  concen¬ 
trate  on 


w  ( X  ,  t )  as  Ae 


lCgi«t 


[(1-i) 


(x-L) 


-fi  (x-L) 
(i-l)e  ‘ 


(4.30) 


We  still  have  A,  which  depends  on  a,  as  an  amplitude  factor  but, 
since  it  multiplies  both  terms,  it  will  play  no  role  in  comparing 
relative  amplitudes  at  different  points  x  in  the  central  region. 
(The  Imaginary  part  of  A  simply  gives  rise  to  a  time  phase  shift.) 
This  is  the  basic  idea  to  be  used  in  identifying  7,  the  damping 
parameter  for  the  beam  itself,  since,  as  we  can  see  from  (4.11), 
is  determined  by  the  coefficients  /> ,  El,  and  7  appearing  in  the 
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equation  (4.10).  So  we  now  consider  Just  the  expression 

,  2(cos  «t  +  1  sin  «t)^sinh  fi^(x-L)  +  1  cosh  (4.31) 

Let  us  set 


S>sin  I  C*  cos  [^“  ®]  =  sin  e  a  c  .  (4.32) 


Then 

if?  -  A 

V  e  ^  ■'  =  !/(£+  iS)  ,  (4.33) 

-  (cos  vSx  *  i  sin  vSx)  .  (4.34) 

Froa  (4.31)and  (4.34)  we  can  see  that 

Re  w(x,t)  =  2  ^cos  wt  cos  i/S(x-L)  sinh  vC(x-L) 

-  sin  «t  sin  i/S(x-L)  cosh  >/C(x-L)  -  cos  «t  sin  vS(x-L)  sinh  i/C(x-L) 

-  sin  (ut  cos  vS(x-L)  cosh  i/C(x-L)j  = 

^cos  «t  cos^vS(x-L)  +  ■jj  cosh  >/C(x-L) 

-  sin  Mt  sin^vS(x-L)  +  sinh  vC(x-L)j  = 

2^^  j^cos^wt  +  fc'S(x-L)  +  +sln  ^«t  -  fc'S(x-L)  -  e~  ^  1 

(4.351 

This  foraula  expresses  the  central  region  motion  in  terms  of  an  "out 
going"  wave  (the  second  term),  attenuated  like  as  it  moves  to 
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the  right,  plus  an  "IncoBlng'*  wave  (the  first  term),  attenuated  like 

as  it  moves  to  the  left.  Prom  the  penultimate  formula  of  (4.35) 
we  see  that,  for  given  x  in  the  central  region,  the  amplitude  of  the 
time  oscillation  is 

j^cos*|^vS(x-L)+-jj  cosh*vC(x-L)  +  s  in*  (x-L  ) s  inh*vC(  x-L  )  j 

2v'^  [cosh*-./C(x-L)  -  sin*[yS(x-L)  +  •  (4.36) 

When  r  (hence  also  C)  »  0  the  last  expression  has  zeros  where 

sin  [vS(x-L)  +  j]  =  ±  1  .  (4.37) 

corresponding  to  "nodes"  of  the  standing  wave  form.  When  7  >  0  there 
are  no  true  nodes;  what  are  left  are  " quas i -nodes "  corresponding  to 
minima  of  (4.36).  The  values  of  these  minima  will  correspond  to  an 

amplitude  factor  times  cosh*>/C(x-L)  -  1 ;  by  comparing  these  minima 
as  they  occur  for  different  values  of  x  in  the  central  region  of  the 
beam  one  may  obtain  an  estimate  for  vC  >  D  for  the  mode  under  study. 
From  (4.11)  we  see  that,  assuming  g,  and  hence  e,  small 


(4.38) 


(4.39) 


(4.40) 


Prom  the  last  formula  we  see  that  a  constant  (i.e.,  froquency 
inctapancfant)  damping  rate,  corresponding  to  7  constant,  is  indicated 
by  the  following  relationship  between  D  and  the  frequency,  u : 
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D  » 


(4.41) 


while  frequBncy  proportional  (sometimes  called  structural)  damping 
corresponds  to  g  constant  so  that 


D 


go 


0( 


(4.42) 


According  to  this  theory,  then,  the  presence  of  structural  damping 
will  be  indicated  by  an  estimated  value  of  D  »  D(<^)  varying  in  pro¬ 
portion  to  as  u  becomes  large. 

In  Figures  4.2  and  4.3  we  show  computed  values  of  the  amplitude 
of  the  time  oscillation  for  two  natural  nodes  of  a  damped  Euler- 
Bernoulll  beam  (these  are  computed  exactly,  rather  than  by  means  of 
the  approximation  (4.36)).  plotted  as  functions  of  the  longitudinal 
bean  parameter.  In  Figures  4.4  and  4.5  we  show  node  shapes  taken  in 
the  laboratory  by  scanning  a  bean  in  forced  vibration  by  means  of  a 
laser  vlbrometer.  The  last  two  figures  do  Indicate  that  the  minima 
corresponding  to  what  we  have  called  the  quasi-nodes  are  smaller  at 
the  right  hand  (free)  end  of  the  beam  than  at  the  left  hand  (clamp 
driven)  end.  The  experimental  data  are  not  of  the  quality  required 
for  reliable  estimation  of  internal  damping  rates  at  the  present  wri¬ 
ting  but  do  indicate  that  improvement  of  our  scanning  techniques  and 
equipment  may  be  expected  to  yield  such  estimates  in  the  near  future. 


C 1  amp 


Figure  4.1  Schematic  for  beam  driven  in  oscillating  clamp 
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